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Abstract 

In this brief note, we consider the variation of the entanglement entropy of a region as the shape 
of the entangling surface is changed. We show that the variation satisfies a Wess-Zumino like 
integrability condition in field theories which can be consistently coupled to gravity. In this case 
the "anomaly" is localized on the entangling surface. The solution of the integrability condition 
should give all the nontrivial finite local terms which can appear in the variation of the entanglement 
entropy. 
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I. INTRODUCTION 



Recently much attention has been paid to the computation of Entanglement Entropy 
in quantum field theories. Entanglement entropy of any quantum system is defined in the 
following way. Given any system we divide it into two subsystems denoted by A and B. If we 
assume that the observables which describe subsystem A and subsystem B, commute among 
themselves then the total Hilbert space of the system can be written as the tensor product of 
the individual Hilbert spaces of subsystems A and B, H = Ha <8) Hb- Let us further assume 
that the total system is described by a density matrix p. We can define a reduced density 
matrix for the subsystem A denoted by pa as pa = Tru B p. The entanglement entropy of 
the subsystem A can be defined as the von Neumann entropy of the reduced density matrix 
Pa, 

S EE {A) = -Tr HA p A lnp A (1.1) 

Entanglement entropy of the subsystem B can also be defined in the same way. Direct 
evaluation of the entanglement entropy in field theories is complicated due to the presence 
of the factor Inp. One can apply the replica trick to compute the trace but in practice 
this can be done only for the simplest field theories or theories with very high degree of 
symmetry. For example, the entanglement entropy is known exactly in 1 + 1 dimensional 
conformal field theories when the subsystems are a line segment and its complement {2,3,6]. 
..t can also be calculated exactly in field theories which have a holographic dual description 
8H15|. For applications to the black hole physics please see 16l-l20|. 

Although the computation of entanglement entropy is extremely difficult, the final answer 
is to some extent universal. For example the leading divergence is in general proportional 
to the area of the entangling surface [l|. In conformal field theories the coefficient of the 
logarithmically divergent piece in the entanglement entropy is also universal, related to the 
conformal anomaly of the theory. There can be similar universal terms in general field 
theories. The finite part of the entanglement entropy can also be universal. So it seems 
that there may be a way of understanding the structure of entanglement entropy which does 
not depend on the details of the field theory under consideration. This is analogous to the 
case of anomalies in field theories. The structure of the anomaly does not depend on the 
details of the field theory. One can determine the structure by solving, for example, the 
Wess-Zumino consistency condition [7]. The coefficients multiplying these anomalous terms 
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certainly depend on the specific field theory. If a similar procedure exists in the case of 
entanglement entropy, then that will be a useful tool for practical computations. 

In this brief note we would like to point out that the entanglement entropy in diffeo- 
morphism invariant field theories 1 have to satisfy quite strong constraints. One can derive 
these constraints in the same way as Ward identities are derived in field theory. The only 
difference is that in this case the path integral is over field configurations which satisfy non- 
trivial boundary condition and as a result the final answer is similar to an " anomalous Ward 
identity" . In particular, it has to satisfy a Wess-Zumino like consistency condition. 

There are some similarities between Wilson loop expectation value and entanglement 
entropy in field theory. When we compute the entanglement entropy using the replica trick 
we essentially compute a path integral over a specific set of field configurations which are 
chosen in such a way that the entangling surface plays a distinguished role. In the case of 
wilson loops things are simpler because a particular loop in the space-time is picked up by 
explicit operator insertion. But in both cases the underlying diffeomorphism invariance of 
the field theory is broken by the choice of the entangling surface or the loop. By this we mean 
the following. The partition function or the quantum effective action of the field theory is 
invariant under the substitution: g — > f*g, where / is the diffeomorphism and f*g is the 
pulled-back metric. But the entanglement entropy or the Wilson loop expectation value is 
not invariant under the substitution, because the diffeomorphism also changes the shape of 
the entangling surface or the loop. So the group of diffeomorphisms will act nontrivially 
on the entanglement entropy or the Wilson loop expectation value. In particular they have 
to transform in such a way that the algebra of diffeomorphism is satisfied and this will 
impose some constraints on the possible forms of the entanglement entropy or Wilson loop 
expectation value. 

II. ENTANGLEMENT ENTROPY 

Let us take the background geometry to be of the form R 1 x M n+1 where R 1 is the time 
direction and M n+1 is the spatial section. We divide M n+1 into two regions, A and B, by 
introducing a codimension one hypersurface £ in M n+1 . We are interested in computing the 

1 By this we mean that the field theory has a conserved stress-energy tensor and so can be consistently 
coupled to gravity. 
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entanglement entropy of the region A. We shall refer to the points belonging to A as points 
inside E. 

Now we apply a spatial diffeomorphism, /, to the background geometry. The background 
metric is of the form: 

ds 2 = g^dafdx" = -dt 2 + j ij (x)dx i dx j (2.1) 

where t is the time coordinate and 7^ is the time independent positive definite spatial metric 
on the constant time slices, M n+1 . Since / is purely spatial it essentially acts only on the 
spatial metric 7^ and its action does not depend on the choice of the time slice. In general 
the diffeomorphism maps the regions A and B to some other regions f{A) and f(B). The 
boundary E between A and B gets mapped to /(E) which is again the boundary between 
f(A) and f(B). We shall assume that the diffeomorphism is continuously connected to the 
identity and hence can be thought of as generated by some vector field. As a result of this 
the diffeomorphism never carries any point of region A or B across the boundary E. In 
short, reflections across the boundary E are not allowed and orientation is preserved by the 
diffeomorphism. This together with the fact that the diffeomorphism is time independent 
has the following interesting consequence. If $ is a field configuration contributing to the 
path integral which computes the entanglement entropy of A, then (f" 1 )*^ is another field 
configuration which contributes to the path integral which computes the entanglement 
entropy of f(A). In fact this leads to the identity SE(f(A),g) = SE{A,f*g). In this 
equation g is the background metric and f*g denotes the pull-back of the background 
metric by the diffeomorphism /. Se{A, g) is the entanglement entropy of the region A with 
the background metric g. If £ is the vector field which generates the diffeomorphism then 
{f*d)fMu = fiV + ^ + Vj,£ M . In our setting the vector field £ is time independent. In the 
following section we shall give a heuristic derivation of the above relation. 



We start with the definition of the entanglement entropy of region A, 

S A = -Tr HA p A lnp A = - lim ^-Tr HA p A n (2.2) 

In practice the trace over the reduced density matrix p A cannot be calculated for arbitrary n 
and so one computes it for integer values of n and then analytically continue it to arbitrary 
values of n. But it should be mentioned that the existence and uniqueness of a proper 
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analytic continuation cannot always be proved. So let us study the quantity Tr BA PA n for 
integer values of n using path integral. 

The matrix element of the reduced density matrix p A is defined as, 

Pa{4>A, <f>' A )\g = Tr H B p\g = J Dg<f> B < 4> A , 4>B\p\4>A, 4>B >g (2.3) 

In the above equation denotes all the fields in the theory. We have decomposed the 
field eigenstate \<j> > on the spatial section M n+1 as \<j> >= \<p A > ®\4>b > and D g (f) B 
is the integration measure over the field configurations (f> B . We have also introduced the 
background metric g on M n+1 and all the integration measures depend on it. Using the 
above definition we can write, 

Tr HA pA 2 \ g = J D g (p A D g( f)' A p A (<pA, <f>' A ) Pa(<I>a> <M ( 2 - 4 ) 

= J D g 4> A D g (f)' A D g (f) B D g (f)' B < <f> A , (f) B \p\(f)' A , (f) B > g < (f)' A , <f>' B \p\<f>A, <f>' B >g 



Now the density matrix of the system is given by, p = where /3 is the inverse tem- 

perature, H is the Hamiltonian and Z is the canonical partition function of the system at 
temperature (3~ l . We can write matrix elements as, 

"<t>(P)=(4>A,4> B ) 
i<P(o)=(4>' a ^b) 

where S is the Euclidean action of the field theory. Now we make a change of variable in 
the path integral, ->■ 4> (= (f^ 1 )*^), 



i r<t>(P)=(4>A,4> B ) 

< <j> A , faW* <pB> g = -= I D g <P e~ s ^ (2.5) 

£ J<h(Q) = (d>',,<bn) 



< MbWaAb > g = t? / D g 4> e- s ^ (2.6) 

z -I m={4>' A ,<t> B ) 



i r (f- i rm=(4>A,<t> B ) 

= - / Dgif' 1 )*^ e- s ^TM 

z ^(/- 1 )*0(O)=(^,<^ B ) 

-, r 4>(f})=u*<t>A,r<t> B ) 

= 1 / Df* g (f) e- s ^g) 



where / is a diffeomorphism which acts only on the spatial section M n+1 and is time inde- 
pendent. In the last line we have used the diffeomorphism invariance of the action and the 
path integral measure. We can write this as an identity, 

< <pA, 4>b\p\4>'a-> 4>B >g = < f*<t>A, f*<pB >f*g (2.7) 
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where we have used the fact that the partition function Z is diffeomorphism invariant. 
Substituting this identity into eqn-(4) we get, 

tt HaPa \ = j D g $ A D g $j) g <\> B D g $ B <r4> A ,r$B\p\r$ A ,r<$>B > f * 9 



(2.8) 

Now we use the diffeomorphism invariance of the measure and write, 



x <r<t> A ,r<t> B \p\f*(f>A,rcf> B >f* 



tt HaPa \ = J D f . g r<i> A D f . g r<i>' A D f . g r<i>BD f . g r<i>' B < > f * 9 

X < />A>/VbIpI/*0A,/>b >f*g 

= Tr Hf _ HA) p f -i {A) 2 \ f * g 

(2.9) 

This equation follows from the facts that f*(p A and f*<pB are field configuration on the 
regions f~ l {A) and f~ l {B) respectively and we are integrating over them with / held fixed. 
It is clear from the above derivation that this is true for any integer n and so it is a reasonable 
assumption that this relation still holds after analytic continuation. So differentiating this 
relation with respect to n we get, 

S E (A,g) = S E (f- 1 (A)J*g) (2.10) 

or 

S E (f(A),g) = S E (A,f*g) (2.11) 

Although we have derived this relation under certain assumptions, our claim is that this will 
hold even in those cases where these assumptions are not valid. In other words, we claim 
that in any diffeomorphism invariant field theory, the entanglement entropy will satisfy the 
relation S E (f(A),g) = S E (A, f*g), where f is any time-independent spatial diffeomorphism 
which is continuously connected to the identity, i. e, generated by some smooth vector field. 
This relation may still be valid for arbitrary diffeomorphism transformations but it certainly 
does not follow from the above logic. 

It also follows from the above equality that the change in the entanglement entropy when 
the shape of the region is deformed is captured completely by the stress-tensor of the theory. 
We showed that changing the shape of the entangling surface is equivalent to changing the 
background metric keeping the shape of the entangling surface fixed. But the response to a 
change in the background metric is captured by the stress tensor of the theory. 
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Now we can expand the entanglement entropy as : 

S E (f(A),g) = S E (AJ*g) = S E (A,g) + J d n+1 x^g8g^^- (2.12) 

If we use the above formula for the change in the metric we get: 

S E (f(A),g) = S E (AJ*g) = S E (A,g) + J d n+1 x^V„(-^^) (2.13) 

It is easy to show by studying diffeomorphisms for which f(A) = A, that the integrand 
has a delta function support on the entangling surface. It is not necessary for / to keep 
A point-wise fixed. In other words the diffeomorphism has to be such that the vector field 
which generates this has to either vanish on the boundary of A or should be tangential to 
the boundary at every point. This is equivalent to the condition f\A) = A. So we can write 

V„rf|^) = M s (2.14) 
y/gog^u 

where n M is the unit outward normal to the surface E and 5s is the delta function supported 
on the surface E 2 .As a result we can write, 

5zS E (A, g) = S E (A, f*g) - S E (A, g) = J F^n, (2.15) 

The right hand side of eqn-(14) could contain derivatives of the delta function of the form 
ft-^V^s, but they can be eliminated by the following argument. If it contains derivatives of 
the delta function then the change in the entanglement entropy under a deformation of the 
entangling surface will contain a term proportional to n M V M £„ evaluated on the surface E. 
Now we can choose £ such that it is vanishing along the surface E but has a nonzero normal 
derivative along E. In that case f(A) = A and using eqn-(13) we get an equation of the 
form, 

J d n xVhF'n^n u V v ^ = Q (2.16) 

This has to be satisfied for all £ which vanish on the surface E. But we can choose the 
normal derivative of £ on the surface arbitrarily even if we keep £ fixed on the surface. So 
the only way this can vanish is if F' = 0. So we do not need any derivative of delta function 



(5s is defined the following integral identity: J M ^/gfS^ = J E VTi/e where M is the space-time, S is a 
codimension one hypersurface, / is a scalar function, g is the background metric ,h is the induced metric 
on the surface S and /s is the scalar function / evaluated on the surface S. 
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on the right hand side. The fact that the R.H.S is proportional to the normal vector, follows 
from the invariance of the entanglement entropy under the reparametrization of the surface 
E. F is some scalar function constructed out of the metric , the normal vector field and 
their derivatives. 

If we are computing the entanglement entropy in a pure state then we must have 
Se(A,q) = Se{B,q). This requires that the function F should be odd under the change 
— > — where is unit normal to the entangling surface. We follow the convention that 
the unit normal points in the direction away from the region for which we are computing 
the entanglement entropy. 

We can see that even eqn-(15) is strong enough to rule out many possible terms in 
the entanglement entropy. For example consider a term in the expression for the entangle- 
ment entropy which can be written as an integral over the entangling surface. If we evaluate 
the left hand side of eqn-(15) on this term then generically this will produce terms containing 
normal derivatives of the vector field on the surface £ and this cannot be integrated out 
because we have a surface integral. If this is the case then that term alone cannot appear 
in the entanglement entropy. So we can either drop that term or if possible add some other 
terms to the entanglement entropy which will cancel the derivatives of the delta function. 
It is easy to show that the Area of the surface £ is a consistent solution of this equation. 
For the surface area contribution, F = K^,, where is the trace of the extrinsic curvature 
of E. We would like to emphasize that we have not proved the area law. We have shown 
that in any diffeomorphism invariant field theory area-law is one of the consistent solutions. 3 

In the next section we shall show that the Wilson loop expectation value satisfies an 
identical type of equation. We choose this as an example because this is somewhat cleaner 
and has many things similar to the entanglement entropy. 

3 This equation is very similar to the loop equation for wilson loops. One can show that the area law in 
the confining phase is a consistent solution of the loop equation, but it is difficult to show that the string 
tension vanishes if there is no confinement, just from the loop equation. 
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III. WILSON LOOP 



In this section we shall show that the Wilson loop expectation value in gauge theory 
satisfies a similar equation. It will make some of our previous arguments clearer. For the 
case of Wilson loops we replace our earlier product geometry with an arbitrary Riemannian 
manifold. For the sake of simplicity we shall also restrict ourselves to Wilson loops in Abelian 
gauge theories. But everything goes through essentially unchanged for non-Abelian loops. 
The expectation value of the Wilson loop is given by, 



were C is a closed contour and S(A, g) is the Euclidean action of the theory. We have 
introduced a background metric g. We further assume that the action S(A, g) and the path 
integral measure D g A are both diffeomorphism invariant. Now we make a change of variable 
in the path integral, 



where / is a diffeomorphism. In the first line we have changed the variable of of integration 
from A to A', in the second line A' is identified with f*A and in the third line we have used 
the diffeomorphism invariance of the action and the path-integral measure. This equation 
can also be written as, 



This equation tells us that changing the shape of the loop is equivalent to changing the 
background metric keeping the loop fixed. The change in the background metric is equivalent 
to stress tensor insertion. In fact it is easy to see that the change is proportional to the 
correlation function of the Wilson loop operator and operators of the form £nV v T Vfl , where 
T^ v is the stress tensor and £ is the generator of the diffeomorphism. The stress tensor 
is covariantly conserved and so the only contribution comes from the contact term which 
is localized on the Wilson loop. The same argument goes through for non-Abelian Wilson 




(3.1) 




(3.2) 



w{f{c),{r l y g ) 



W(f(C),g) = W(C,f*g) 



(3.3) 
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loops. For simplicity let us consider only non self-intersecting or simple loops. Since we 
are only considering diffeomorphisms, simple loops get mapped to simple loops and the 
orientation is preserved by the diffeomorphism according to our assumption. 
We now expand the Wilson loop expectation value as, 

W(f(C),g) = W(C, f*g) = W(C, g) + f d n x^g £ h V a (^L^) (3.4) 

where £ is the vector field generating the diffeomorphism /. So as long as the diffeomorphism 
does not move the loop, i.e, }\C) = C ,we get 

-2 SW 

d n x v ^aV a ( — — ) = (3.5) 

y/gdgab 



then by standard arguments, we get 

_ ,-2 5W 



a Vd 6 9< 



ab 



F^n b (l) S c (3.6) 



where the n^-s are (n — 1) orthonormal vectors all of which are orthogonal to the loop 
C. The right hand side is proportional to the vectors rii because the expectation value is 
invariant under the reparametrizations of the loop. F^-s are (n — 1) scalar functions built 
out of the metric, the normal vectors and their derivatives. F l must transform in the vector 
representation of the local SO(n — 1) group which acts on the set of orthonormal vectors, 
rii. This follows from the redundancy in our choice of the normal vectors which are defined 
only up to local SO(n — 1) transformation. The change of W(C, g) under a deformation of 
the loop C can be written as, 

W(f(C),g) = W(C,f*g) = W(C,g) + [ (3.7) 

Jc 

A. W Vs. InW 

So far we have not talked about the nature of the functions F l . In general they will not be 
local, i.e, a scalar function built out of the metric, the normal vectors rii and their derivatives 
evaluated at a given point of the space. In particular the variation of the Wilson loop will 
produce functions which will not have this property in general. But InW is expected to be 
much better behaved in this respect. Our previous analysis remains unchanged if we replace 
W with InW and so we can write, 

lnW(f(C),g) = lnW(CJ*g) = lnW(C,g) + / (3.8) 

Jc 

10 



where the functions F l are expected to contain local terms. It can also contain non local 
terms whose existence cannot be ruled out by this procedure. 

IV. WESS-ZUMINO LIKE INTEGRABILITY CONDITION 

In the previous sections we have studied entanglement entropy and Wilson loop expec- 
tation values in arbitrary field theories. We have shown that they both behave similarly 
when the shape of the entangling surface or the loop is deformed. In both cases we have 
written down the variation in terms of some unknown functions. But is there any way to 
determine them? Entanglement entropy (Wilson loop expectation value) or more precisely 
the function F(F % ) has to satisfy an integrability condition. Eqn-(2) gives the change of 
entanglement entropy under a diffeomorphism. But given two regions there exist in general 
more than one (infinitely many) diffeomorphisms connecting them and all of them should 
give the same change in entanglement entropy. This is equivalent to the condition that the 
algebra of diffeomorphism has to be obeyed. It leads to the following condition, 

The vector fields in the commutator appear in the reverse order because we are considering 
active diffeomorphism. One can write eqn(5) in terms of the function F as 

F) - F) = Q([6, F) (4.2) 

where Q is the integral 

Q(£,F) = Jd^Vhn^F (4.3) 

Here h is the induced metric on the entangling surface S. This is our main result. We have 
checked that F = K^, which comes from the area term in the entanglement entropy satisfies 
this integrability condition. It is easy to check that any arbitrary function F will not satisfy 
this integrability condition. So this condition is nontrivial. Similarly we can write down the 
same consistency condition for the Wilson loop. 

A. Comments on The Solution of The Integrability Condition 

We hope that the integrability condition we have stated above will be solvable in the 
same way as the Wess-Zumino consistency condition. But the difference with the standard 
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anomaly answer is clearly visible here. In our case the "anomaly" is localized on a finite 
codimension (entangling) surface. Also in the case of anomaly in field theory, the gauge 
variation of the quantum effective action is always local and finite. But in this case it is not 
clear if the function F(F l ) is always local. In general it will contain both local and non-local 
terms. Also the entanglement entropy and the Wilson loop expectation values are both 
ultraviolet divergent quantities. If we assume that the field theory has been regularized in 
such a way that diffeomorphism invariance is respected, then this should not cause any 
problem. The interesting point is that one can determine the structure of the local terms 
in F(F l ) by solving the consistency condition. There will be two kinds of solutions. The 
trivial solutions will be those which can be obtained by applying a diffeomorphism to the 
integral of a local function. The integral can be a surface integral on E or a volume integral 
defined on the region bounded by S. But the interesting solutions are those which cannot 
be obtained in this way and presumably there is only a finite number of solutions with 
this property 4 . In other words, the nontrivial solutions of the integrability condition define 
the nontrivial local terms which can appear in the variation of the entanglement entropy. 
In many cases the variation of the entanglement entropy is more well defined than the 
entanglement entropy itself because the ultraviolet cutoff may cancel in the variation. The 
solution of the integrability condition we have stated above will give the nontrivial local 
terms which can appear in the variation. It is important to note that in order to solve this 
integrability condition one does not need to know about the field theory, except that the 
specific numerical coefficients multiplying the purely geometric expressions will depend on 
the theory. It may happen that some of the local contributions in F are actually protected. 
This may provide us with some examples of protected quantities in non-supersymmetric 
field theory. The other point is that this method does not require introduction of a 
background geometry with conical defect and so the calculation is under control at every 
stage. 



There are only a finite number of them because in general the solutions of this problem are the cohomology 
classes of a BRST like operator. 
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